PROCEEDINGS 

OF SCIENCE 



Reconstructing Dark Matter Properties via 
Gam ma- Rays with Fermi-LAJ 



Nicolas BERNAL 

Bethe Center for Theoretical Physics and Physikalisches Institut, Universitat Bonn 
Nussallee 12, D-53115 Bonn, Germany 

Centra de Fisica Tedrica de Partwulas, Instituto Superior Tecnico 
Avenida Rovisco Pais, 1049-001 Lisboa, Portugal 



E-mail: nicolas@th. physik . uni-bonn . de 



We study the capabilities of the Fermi-LACT instrument for identifying particle Dark Matter pro- 
perties as mass, annihilation cross section and annihilation channels, with gamma-ray observa- 
tions from the Galactic Center. For the potential Dark Matter signal, besides the prompt gamma- 
ray flux produced in Dark Matter annihilations, we also take into account the flux produced by 
inverse Compton scattering of the electrons and positrons generated in Dark Matter annihilations 
off the interstellar photon background. We show that the addition of this contribution is crucial in 
the case of annihilations into e + e~ and pairs. In addition to the diffuse galactic and extra- 

galactic background, we also consider the full catalog of high-energy gamma-ray point sources 
detected by Fermi. The impact of the degeneracies between the different Dark Matter annihila- 
tion channels has been studied. We find that for Dark Matter masses below ~ 200 GeV and for 
typical thermal annihilation cross sections, it will be possible to obtain stringent bounds on the 
Dark Matter properties. 
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1. Introduction 

If Dark Matter (DM) is detected and identified, the measurement of its properties like mass, 
annihilation cross-section and annihilation channels plays a central role in the determination of 
the particle nature of the DM. It will allow us to constrain models of particle physics beyond 
the Standard Model, for instance supersymmetry and universal extra dimensions. Furthermore 
a convincing DM discovery may require consistent signals in multiple experiments in multiple 
channels (direct, indirect, collider). We discuss the capabilities of the Fermi-LAT instrument for 
identifying particle DM properties with gamma-ray observations from the Galactic Center (GC). 

The differential intensity of the photon signal from a given observational region in the galactic 
halo from the annihilation of DM particles has different possible origins: internal bremsstrahlung 
and secondary photons (prompt) as well as Inverse Compton Scattering (ICS). External brems- 
strahlung and synchrotron emission also contribute to the photon flux; however, for the energies 
of interest here and for typical DM masses, both bremsstrahlung and synchrotron emission are 
expected to be subdominant with respect to ICS. For the sake of simplicity we will neglect these 
sources in what follows. 

The differential flux of prompt gamma-rays from DM annihilations and coming from a direc- 
tion within a solid angle AO. is given by 

— * (E r ,AQ. =^4£-lBR ; -— / da I p{r(s,a)) 2 ds, (1.1) 
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where (crv) is the total thermally averaged annihilation cross section, m x the mass of the DM 
particle, BR, the annihilation fraction into channel i, dN'y/dEy the differential gamma-ray yield of 
standard model particles into photons of energy Ey, p (r) the DM density profile and r the distance 
from the GC. Here we will focus on the NFW halo profile [P; the dependence on the DM halo 
profile has been studied in reference 

An abundant population of energetic electrons and positrons produced in DM annihilations 
either directly or indirectly from the hadronization, fragmentation, and subsequent decay of the 
SM particles in the final states, gives rise to secondary photons at various wavelengths via ICS off 
the diffuse radiation fields in the galaxy. We approximate this photon background as a superposition 
of three black-body spectra consisting of the CMB, the optical starlight and the infrared radiation 
due to rescattering of starlight by dust The differential flux of high energy photons produced 
by the ICS processes is given by [Q] 

(d<f> \ 1 1 t f r m x dn 
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where &(Ey, E) is the differential power emitted into scattered photons of energy Ey by an electron 
with energy E. The minimal and maximal energies of the electrons are determined by the electron 
mass m e and the DM particle mass. The quantity dn e /dE is the electron plus positron spectrum 
after propagation in the Galaxy, which will differ from the energy spectrum produced at the source. 
We determine the propagated spectrum by solving the diffusion-loss equation that describes the 
evolution of the energy distribution for electrons and positrons assuming steady state [Qj. Regarding 
the propagation parameters (like diffusion coefficient, energy losses and thickness of the diffusion 
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zone), we take their values from the commonly used MED model Again, the dependence on 
the propagation model has been studied in reference 

There are three main components contributing to the high-energy gamma-ray background: the 
diffuse galactic emission has been estimated by taking the conventional model of the GALPROP 
code On the other hand, another source of background particularly important when looking 
at the GC is that of resolved point sources. We consider all the point sources detected by the first 
1 1 months of Fermi-LAI [f7|] lying in the region of interest. Finally, for the isotropic extragalactic 
gamma-ray background we used the recent measurements by the Fermi-LAI collaboration 
For a 10° x 10° region around the GC, the diffuse galactic emission dominates below ~ 20 GeV. 
Above that value, the emission coming from point sources is the most important. The isotropic 
extragalactic gamma-ray background is at the percent level. 

The Large Area Telescope (Fermi-LAI) is the primary instrument on board of the Fermi 
Gamma-ray Space Telescope. It performs an all-sky survey, covering a large energy range for 
gamma-rays, with an effective area ~ 8000 cm 2 and a field of view of 2.4 sr. In the following ana- 
lysis, we consider a 5-year mission run, and an energy range from 1 GeV extending up to 300 GeV. 
We divide this energy interval into 20 evenly spaced logarithmic bins. In order to maximize the 
signal-to-noise ratio, it has been pointed out that for a NFW profile the best strategy is to focus on 
a region around the GC of ~ 10° x 10° [g]. Hence, this is our choice. 



2. Reconstructing Dark Matter properties 

Once gamma-rays are identified as having been produced in DM annihilations, the next step 
concerns the possibilities of constraining DM properties [[Tc|]. In figure [j] we depict the Fermi- 
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Figure 1: Fermi-LAI abilities to constrain DM properties. We consider DM annihilation into a pure bb 
final state and two DM masses: m x = 80 GeV (left panels) and m x = 270 GeV (right panels). Dark blue 
(light orange) regions represent 68% CL (90% CL) contours. We assume a 10° x 10° observational region 
around the GC, a NFW DM halo profile, the MED propagation model and (ctv) = 3 • 10~ 26 cm 3 s -1 . The 
black crosses indicate the values of the parameters for the simulated observed "data". 
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LAT reconstruction prospects after 5 years for DM annihilation into a pure bb final state re- 
constructed as either t + t~ or bb and two possible DM masses: m x = 80 GeV (left panels) and 
m x = 270 GeV (right panels). We also assume DM particle with a typical thermal annihilation 
cross section (av) = 3 • 10~ 26 cm 3 s -1 , the MED propagation model, a NFW DM halo profile 
and a 10° x 10° observational region around the GC. These benchmark points are represented in 
the figure by black crosses. The dark blue regions and the light orange regions correspond to the 
68% CL and 90% CL contours respectively. The different panels show the results for the planes 
(m x , (av)), (BR T (/,), (av)) and (m x ,BR T ^)), marginalizing with respect to the other parameter 
in each case. BR T ( fo ) = 100% (0%) corresponds to an annihilation into a pure T + T~ (bb) final 
state. For the first model chosen in figure [j] (left panels), m x = 80 GeV and in general for light 
DM masses, the reconstruction prospects seem to be promising, allowing the determination of the 
mass, the annihilation cross section and the annihilation channel at the level of ~20% or better. On 
the other hand, for heavier DM particles, the regions allowed by data grow considerably worsening 
the abilities of the experiment to reconstruct DM properties. This is shown for the second model in 
figure [I] (right panels), m % = 270 GeV. In this case, Fermi-LAT would only be able to constrain the 
DM mass to be in the range ~ (30 — 500) GeV and determine the annihilation cross section within 
an order of magnitude. Let us note the appearance of a second spurious minima corresponding to 
a lighter mass m x ~ 60 GeV, an annihilation cross section (av) ~ 10~ 26 cm 3 s _1 and annihilating 
mainly (~ 75%) on t + t~ pairs. 
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Figure 2: Fermi-LAT abilities to constrain DM properties. We assume the measured signal is due to DM 
annihilating into but the fit is obtained assuming DM annihilates into either T + T~ or bb. We assume 

two DM masses: m x = 50 GeV (left panels) and m x = 105 GeV (right panels). Dark blue (light orange) 
regions represent the 68% CL (90% CL) contours. See the text for the rest of the parameters. The black cross 
in the left-top panel in each plot indicates the values of the parameters for the simulated observed "data". 
Note that the other panels have no cross as they lie outside the parameter space of the simulated observed 
"data". The squares indicate the best-fit point. 

In figure ^| we assume that DM actually annihilates into jU + /i~ pairs, but we analyze the data 
assuming DM annihilations into either t + t~ or bb, for (av) = 3 • 10~ 26 cm 3 s _1 and for two DM 
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masses: m x = 50 GeV (left panels) and m x = 105 GeV (right panels). Again, the black crosses 
indicate the values of the parameters for the simulated observed "data"; the squares indicate the 
best-fit points. Naively, one would expect that the jU + jU~ (leptonic) channel is identified as being 
closer to the t + t~ (leptonic) channel than to the bb (hadronic) channel. Let us remember that 
DM annihilation channels are commonly classified into two broad classes: hadronic and leptonic 
channels. Leptonic channels typically give rise to a harder spectrum and, in particular for e + e~ and 
, the cutoff is be very sharp and around a maximum energy (i.e. the mass of the DM particle 
in the case on annihilating DM). Contrary to what was expected, the reconstructed composition of 
the annihilation channels tends to be dominated by bb, instead of t + t~. Indeed, the contribution 
due to ICS in the case of the (and also the e + e~) channel could substantially alter the 

different prompt spectra. Hence, when taking into account the contribution of ICS to the gamma- 
ray spectrum, the annihilation channels cannot be generically classified as hadronic or leptonic, as 
DM annihilations into jU + jU~ pairs are better reproduced with the bb channels than with the t + t~ 
channel. 
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Figure 3: Fermi-LNT abilities to constrain DM properties. We assume the measured signal is due to DM 
annihilating into and the fit is obtained assuming DM annihilates into or bb. We assume 

ICS+prompt photons (left panels) or only prompt photons (right panels) for the reconstructed signal, for 
m x = 50 GeV. Dark blue (light orange) regions represent the 68% CL (90% CL) contours. See the text for 
the rest of the parameters. The black crosses indicate the values of the parameters for the simulated observed 
"data". The squares in the right panels indicate the best-fit point. 

The results just discussed can be illustrated in a different way by analyzing the simulated 
observed signal "data" from DM annihilation into pL + pL~ assuming DM annihilates into either 
or bb. This is depicted in figure || where we show the results for the case that we try 
to reconstruct the full signal (prompt and ICS) generated by a 50 GeV DM particle adding the 
ICS contribution (left panels) or with only prompt photons (right panels). As can be seen in the 
left panels, if ICS is taken into account, DM properties can be reconstructed with good precision. 
However, if the ICS contribution is not added to the simulated signal events (the simulated observed 
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"data" always has the ICS included), DM annihilation into a pure channel would be excluded 

at more 90% CL, providing thus a completely wrong result. 

3. Conclusions 

In this work we have studied the abilities of the Fermi-LAT instrument to constrain Dark Mat- 
ter properties by using the current and future observations of gamma-rays from the Galactic Center 
produced by DM annihilations. Unlike previous works, we also take into account the contribution 
to the gamma-ray spectrum from ICS of electrons and positrons produced in DM annihilations off 
the ambient photon background. We show that the inclusion of the ICS contribution for hadronic 
channels and for the T + T~ channel does not give rise to important differences in the reconstruction 
process. This is not the case if DM annihilates into the jJ. + jJ.^ and the e + e~ channel. In this latter 
case, adding the ICS contribution to the prompt gamma-ray spectrum turns out to be crucial in 
order not to obtain completely wrong results. 

On the other hand, we found that for Dark Matter masses below ~ 200 GeV and for typical 
thermal annihilation cross sections, it will be possible to obtain stringent bounds on the Dark Matter 
properties such as its mass, annihilation cross section and annihilation channels. 
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